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DIFFUSION COMPRESSOR BLADING 
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# 
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ABSTRACT 

A method Is presented for automating compressor blade design using 
numerical optimization, and Is applied tc the design of a controlled diffu- 
sion stato*’ blade row. A general purpose optimization procedure is em- 
ployed, which Is based on conjugate directions for locally unconstrained 
problems and on feasible directions for locally constrained problems. 

Cojpled to the optimizer Is an analysis package consisting of three analysis 
programs which calculate blade geometry, Invlscld flow, and Made surface 
bounoary layers. 

The optimization concepts are briefly discussed. Selection of design 
objective and constraints is described. The procedure for automating the 
design of a two-dimensional blade section Is discussed, and design results 
are presented. 
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NOMENCLATURE 


%1 

AC1, 601,001,001, 

polynomial coefficients for blade 


AC2,BC2,CC2,0C2 

angle distribution expression 

%le 

ATl.BTl.CTl.DTl 

polynomial coefficients for blade 

AT2,BT2,CT2,0T2 

thickness distribution expression 


C 

chord 

s m2 


objective function 

constraint functions 


Hi 

incompressible form factor 

s m2e 

Hi,crit 

critical value of incompressible 
form factor at which turbulent 


boundary layer separates, 

H i crit * this study 

T 

L.E. 

leading edge of blade 


KCS 

angle with respect to meridional 



direction of blade mean line mid- 
way between transition location 

t l.t2 


and trailing edge 

TMX 

KICR 

angle with respect to meridional 

T.E. 


direction of blade mean line at 

U 


leading edge 

V 

K0CR 

angle with respect to meridional 

X 


direction of blade mean line at 
trailing edge 

XCH0R0 

KTC 

angle with respect to meridional 
direction of blade mean line at 

ZM 


transition location 

a* 

X 

search direction 

K 

si 

distance from transition location 


along mean line, front segment 

Subscripts 

s 2 

distance from transition location 

along mean line, rear segment 

1 

0 

^Member ASME. 


c 

m 


n 

ps 

ss 

t 


distance from maximum thickness 
location along mean linf, front 
segment 

distance from maximum thickness 
location to blade leading edge 
along mean line 

distance from maximum thickness 
location along mean line, rear 
segment 

distance from maximum thickness 
location to blade trailing edge 
along mean line 
distance from leading edge to 
intersection of two polynomial 
segments describing mean in el 
chord 

thickness of blade, front and roar 
segments respectively 
maximum thickness/chord 
trailing edge 
stream function 
surface velocity 
vector of design variables 
meridional projection of blade 
chord 

distance from leading edge to 
maximum thickness location/chord 
move parameter 
blade angle 


front segment 

rear segment 

number of constraints 

number of design variables 

pressure surface 

suction surface 

transition 


INTRODUCTION 

Throughout the history of compressor technology * 
blade shapes have been specified by geometric fam- 
ilies or classes. For the most part, these families 
have been derived from early wing shapes and improved 
bv empiricism, or have been directly specified from 
simple geometric shapes such as circular arcs and 
parabolas. 

During the past decade, computational methods 
for the calculation of flow through compressor blade 
rows have advanced substantially, as have computer 
speeds. With these advances has come the capability 
to rapidly design and analyze flow over arbitrary 
blade shapes. Indeed, at the present time, thesu 
analysis methods a-.i being synthesized into computer- 
aided design systems. In most cases these systems 
are "manual", 1,c., non-automated. Because of the 
great flexibility in choice of blade shape, the 
design process can become quite cumbersome &nd repe- 
titive unless automated in some fashion. One of the 
most attractive methods for automating the design 
process is numerical optimization. Much progress has 
been made recently in bringing the technique to bear 
on engineering problems, particularly in the field of 
Aeronautics (1). Of the many numerical optimization 
algorithms in existence, the one used in iUf. 1 and 
described in Ref. 2 with its control program (3) is 
sufficiently general and user-oriented to be of par- 
ticular interest. It Is used in the work reported 
herein, and is coupled to analysis programs which 
calculate blade shape, the inviscid flow field, and 
the boundary layer for a two-dimensional blade 
section. 

With the advent of arbitrary blade shapes, the 
concept of controlling velocity diffusion (and con- 
sequently boundary layer growth) on tlve suction sur- 
face has received increasing attention. In the 
transonic flew regime, such blading has generally 
been referred to as "supercritical’ blading" since the 
local supersonic flow is controlled as well as the 
boundary layer growth, in the subsonic regime the 
blading is often simply referred to as ''controlled 
diffusion". Methods of analysis have generally been 
inverse, in which a velocity distribution of a gen- 
eral Stratford type (a) is prescribed at the outset, 
and a blade shape derived from it (5 and 6). 

The problem addressed in the present work is the 
redesign of a high-subsonic stator blade row uti- 
lizing a controlled diffusion blaae shape. The ana- 
lytical methods are direct rather than Inverse. 

A blade shape is initially prescribed and aerodynamic 
performance calculated. Perturbations on the blade 
shape are effected and aerodynamic performance re- 
calculated until specified conditions are met. The 
resulting velocity distributions over the suction 
surface of the blade are also of the general 
Stratford type, but in this case are controlled by 
constraints imposed on the geometric and aerodynamic 
parameters. 

The subject stator row uses the same flow path 
and velocity triangles as the first stage stator of 
the NASA Two-Stage Fan (7). The original design was 
highly successful, showing a first stage peak adia- 
batic efficiency of 87.0 percent, and a remarkably 
low radial distribution of loss across the stator. 
Consequently, significant improvement in performance 
with controlled diffusion blading cannot be expected, 
nor is that the purpose of the present work. The 
principal objective of the work presented herein is 
to develop and demonstrate the feasibility of an 
automated design procedure based on numerical opti- 
mization, Experimental evaluation of the resulting 


design is planned for both a single-stage environment 
and a two-dimensional cascade (midspan blade section). 

ANALYSIS METHODS 

Blade Section Geometry 

"The blade section geometry is generated by a 
polynomial element program. This program has been 
extracted from the NASA Design Program which Is a 
streamline curvature design procedure (0). Blade 
section nomenclature is presented In Fig. 1, The 
meanline of the blade is described by two polynomial 
segments, each of which can be specified by up to a 
quartic polynomial. The polynomial is a fit of local 
mean-line blade angle against mean-line distance. 

The fraction of chord from the leading edge at which 
the two polynomial segments join is referred to as 
the transition location, T, The polynomials may be 
fitted beginning from the transition location and 
fitting toward the leading edge and trailing edge 
respectively, or they may be fitted beginning at the 
leading and trailing edges and fitting toward the 
transition location. In this report, fitting from 
transition location toward leading and trailing edges 
for each segment is the mode of operation. Note 
that si and s? are both positive in this 
mode. 

The expression for blade angle distribution is 
given by Eq. (1): front segments 

tcj « K t + AC1 x s 1 + BC1 x + CC1 x s l + DC1 x sj 


(1A) 


Rear segments 


Kg « + AC2 X Sg + BC2 x s| + CC2 x Sg + DC2 x Sg 


(IB) 


A typical distribution of blade angle is shown 
in Fig. 2(a). The blade angle at the transition 
location is designated as. KTC. When fitting from 
the transition location, the fit is effected from 
KTC to the blade angle at inlet, K1CR (front seg- 
ment), and from KTC to the blade angle at outlet, 
KOCR (rear segment). Noted on the figure is a 
parameter calculated Internally by the program, K£S, 
which is the blade angle midway between the transi- 
tion location and the trailing edge. This parameter 
will be of importance in later discussion. 

The distribution of blade thickness about the 
mean-line is also specified by two polynomials, both 
of which may be quartics. The thickness is added 
symmetrically on either side of the meanline. The 
fit is made from the maximum thickness location 
toward the leading and trailing edges for front and 
rear segments respectively. The leading edge and 
trailing edge may be specified as circles or 
ellipses. Circles only were used in this design. 

The equations for thickness distribution are given 
as Eq. (2). Instead of being linear, the first term 
is of a square root form, wit lew enables simulation 
of 65-series blades, if desired. 
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Front segment,* 

*•1 " TT * AT1 x ^*mle ~ s ml " V^mle * 2 ^ 

- BT1 x sj a - CT1 x 4 - 0T1 x ^ (2A) 


Rear segment,* 


. BT2 X 4 - CT2 X s 3 , - 0T2 X sj 2 (2B) 


Potential Flow Solut i o n 

The potential flow about the blade section In 
the two-dimensional, blacie-to-blade plane Is cal- 
culated by the method developed by Katsanis, 

TSQNIC (9). The program solves the stream func- 
tion equation by finite difference techniques for 
the subsonic, compressible flow regime. It Is 
necessary to specify as Input the fluid proper- 
ties, inlet total temperature and density, weight 
flow, blade geometry, Inlet and outlet flow 
angles, finite difference mesh, and a meridional 
distribution of streamtube height and total pres- 
sure loss. In the design presented herein, a 
linear distribution of streamtube height and esti- 
mated total loss was utilized. 

Because the nature of the equations dictates 
that the solution be of a boundary value type, the 
outlet flow angle must be specified on the down- 
stream boundary. This efrectively sets the Kutta 
condition. Since this condition is related to one 
of the constraints chosen for the optimization 
process, its discussion will be reserved until 
later. 

Boundary Layer Calculations 

Blade surface boundary 1 layers were calculated 
using the program developed by McNally (10). in 
addition to the surface velocities, required input 
includes upstream flow conditions, fluid proper- 
ties, and blade surface geometry. Among the out- 
put provided by the program are the conventional 
boundary layer thicknesses, form factors, wall 
friction coefficient, and momentum thickness 
Reynolds number. 

The program uses integral methods to solve 
the two-dimensional compressible laminar and tur- 
bulent boundary layer equations in an arbitrary 
pressure gradient. Cohen and Reshotko's method 
(11) is used for the laminar boundary layer, 
transition is predicted by the Schlictlng-Ulrich- 
Granville method (12), and Sasman and Cresci's 
method (13) is used for the turbulent boundary 
layer. 

A boundary layer which is initially laminar 
may proceed through normal transition to a turbu- 
lent boundary layer, or it may undergo sou farm 
of laminar separation before becoming turbulent. 

To provide flexibility for analyzing this be- 
havior, several program options are available to 
the user. The calculations may proceed from a 
laminar boundary layer through transition to tur- 


bulent calculations. However, if laiiiihar separa- 
tion is predicted before transition, the turbulent 
calculations may be started by specifying a factor 
by which the last calculated value of momentum 
thickness is multiplied (this value is commonly 
chosen to be 1.0 to satisfy conservation of 
momentum). This new momentum thickness and a 
value for form factor based on the last calculated 
momentum thickness Reynolds number are used as 
initial values for the turbulent calculations, 

Optimization Program 

The optimization algorithm in Fortran code 
Is known as C0NM1N, and Is reported In Ref. 2. 

A general purpose control program known as COPES 
is coupled to the algorithm (3). 

The general mathematical representation of a 
numerical optimization problem is stated as: 


Minimize OBJ w F(5T) 

subject to Gj(5C) < 0. 

x) < x 1 < X 

J is a vector consisting of the design variables. 

Xf and XV are the lower and upper bounds 
on the design variables and are referred to 
as side constraints. OBJ is the objective func- 
tion. If the designer wishes to maximize a func- 
tion, OBJ may be defined as the negative of the 
function. Gj{¥) set the constraint functions 
which the design must satisfy. When Gj(X) < 0, 
it Is said to be inactive; when Gj(X') > 0, it is 
violated. When it is within a tolerance band about 
zero, it is active. F0f) and Gj(lf) may be im- 
plicit or explicit functions of the design vari- 
ables K, but must be continuous. (Note: this 
should be carefully considered when formulating 
these functions when they are calculated from finite 
difference solutions or at discrete stations.) 

An initial design vecto**, IT, is specified by 
the user. It may be feasible or Infeasible, i.e., 
if it satisfies the inequalities of Eq. (3), it is 
feasible. If a feasible Initial design can be 
found, it is usually more efficient to begin with 
it, at least for the types of problems discussed 
herein. An Iteration process is then begun which 
follows the recursive relationship: 


TP * 1 • (A) 

q is the iteration number; the vector is the 
search direction In the n-dimensional space; and 
the scalar a* (move parameter) defines the dis- 
tance of travel in direction "S, and is found by 
interpolation, 

The search direction "S is initially obtained 
by moving in the direction of steepest descent 
(negative gradient of the objective function) with- 
out violating constraints. The procedure is then 
repeated using a conjugate direction algorithm in 
determining a new search direction. Whenever a 
constraint is encountered, a new search direction is 
found using Zoutendijk's Methoo of Feasible Direc- 
tions. An optimum has been achieved when no search 
direction can be found which will further reduce the 
objective function without violating a constraint. 


l,m 
1 , n 


(3) 
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OPTIMIZATION OF STATOR BLADE SECTION 

Formulation of a specific optimization problem 
Involves choice of an objective function (the quan- 
tity to be optimized), choice of constraints, and 
choice of deslqn variables. In the present design 
problem, optimization of a two dimensional stator 
blade section was performed at the 90 percent span 
from tip section. This location represented the 
most difficult design problem as measured by blade 
loadlna requirements, 

Results of preliminary calculations of an 
initial blade shape which meets the specified 
velocity triangles at the 90 percent span location 
are shown In Figs. 2 to 4. Figure 2 shows the dis- 
tribution of blade angle and blade thickness along 
the mean-line, and the corresponding blade shape, 
Figure 3 is the surface velocity distribution 
corresponding to the blade shown In Fig. 2. And 
Fig. 4 represents the incompressible form factor 
distribution along the suction surface obtained from 
the boundary layer calculations. 

The Initial blade design was essentially an 
arbitrary choice. The blade angle distribution and 
thickness distribution plots were determined by 
running the blade geometry program In a graphics 
mode. In this mode, the distributions can be 
generated by curvefitting through points which are 
input by the user. The process Is therefore Intui- 
tive, and guided by experience. The only restric- 
tion to the process is the desirability, with regard 
to optimization theory, that the design be feasi- 
ble. Note that for the initial design selected, the 
turbulent boundary layer separates at 64 percent of 
chord. 

A properly designed controlled diffusion blade 
should experience no suction surface boundary layer 
separation. This criterion is incorporated 1,ne the 
objective function. The following penalty function 
type of objective function proved to be the r os'- 
successful , 

OBJ - FORMAX - XSEPOX (5) 

FORNAX is the maximum incompressible form factor 
(Hi) occuring 'ver the rear portion of the blade, 
and XSEPOX is the separation location of the turbu- 
lent boundary layer expressed as a proportion of 
chord length. 

OBJ was minimized. Reducing FORMAX acts to 
increase the separation location, XSEPOX. Simulta- 
neously reducing FORMAX and increasing XSEPOX acts 
to reduce OBJ. 

Design Variables 

'Nine design variables have been selected, all 
of which describe the geometry of the blade. These 
variables are: 

T, transition location of two mean-line 
polynomials; 

ZM, maximum thickness location; 

KOCR, the exit blade mean-line angle (deg.); 

AC1.BC1, first two coefficients of front seg- 
ment mean-line polynomial; 

AC2,BC2,CC2,DC2, all four coefficients of rear 
segment mean-line polynomial. 

The velocity triangles are fixed for the blado 
section, thus f'vxing the loading or overall velocity 
diffusion across the blade row. By allowing KOCR 
to vary, the blade camber angle is allowed to 
change. The trailing edge or Kutta-type condition 
is controlled through a constraint described below. 


Incidence angle Is not allowed to vary, It is 
fixed at the value used in the original design (7) 
Simply as a designer's preference. If It should 
become desirable or necessary to allow incidence 
angle to vary, It can easily be Incorporated by 
including KICR as a design variable, and retaining 
the same velocity triangle Information. 

Although the maximum thickness location, ZM, is 
allowed to vary, the coefficients of the blade 
thickness polynomial are held fixed at the values 
used in the preliminary (initial) blade shape, 
strictly as a designer's choice. 

Each of the above variables is allowed to vary 
within user-selected limits. The upper and lower 
bounds for each are listed below, and in optimiza- 
tion theory are referred to as side constraints. 

The side constraint values are chosen as a result of 
experience. Beyond certain values, a realistic 
blade shape will not result, In addition, certain 
extreme combinations of variables may cause 
convergence problems in the geometry program, and so 
are best avoided. 



Lower 

Upper 


Bound 

Bound 

Transition Location/Chord, T 

.20 

.40 

Max Thickness Loc, /Chord, ZM 

,35 

,55 

Outlet Blade Angle, KOCR 

-10. 

-2. 

A1 1 Coefficients 

1 . E— 15 

l.E+15 


Constraints 

Fi«e constraint functions are specified, all 
being implicit functions of the design variables. 

Two constraints are variables calculated Internally 
to the geometry program and control the blade angle 
distribution. By controlling the blade angle dis- 
tribution, a controlled diffusion type shape to the 
surface velocity distribution can be Insured. The 
constraints are represented in Fig. 2(a), and are 
KTC and KCS, which were previously described. 

They were allowed to vary between the following 
bounds: 

Lower Bound Upper Bound 

KTC 32,0 46.08 

KCS -4.0 11.05 

The remaining constraints arc calculated in the 
inviscid flow program (9) and are described with 
reference to a surface velocity distribution. Since 
TSONIC is principally a subsonic calculation proce- 
dure, the maximum surface velocity on the suction 
surface is constrained to the subsonic flow regime, 

A search procedure locates the maximum suction sur- 
face velocity, This is nondimensionalized by the 
inlet freestream velocity and the ratio is defined 
as the constraint. The upper bound was set to be 
equivalent to Mach 1 condition. The lower bourn# is 
set equal to an arbitrary, small number. :• 

Because no boundary layer calculations are made 
on the pressure surface of the blade, a constraint 
is applied to control the velocity diffusion on the 
pressure surface. Preliminary calculations were 
made of typical blade shapes, and a pressure surface 
velocity diffusion (V,max/V,min) of 1.65 was deemed 
to be a sufficiently safe upper bound, Subsequent 
to the optimization calculation, the pressure sur- 
face boundary layer of the optimized blade is also 
calculated to verify that it was truly free from 
separation. The lower bound of this constraint is 
set equal to an arbitrary, small number. 
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Him final constraint Is chosen to set the 
trailing edge condition, the condition equivalent to 
a Kutla condition, this Is also equivalent to 
setting deviation angle. In the present work, out- 
le v flow angle is fixed, so whatever value is taken 
on by blade outlet ang^e, KOCR, sets deviation 
angle. Experience with some of the conventional 
families of blades, supplemented with detailed 
analyses, permitted guidelines to be set for esti- 
mating deviation for those blade families (14). But 
such experience Is lacking for controlled diffusion 
blading, which Is arbitrary In shape. For conven- 
tional blading, setting the deviation angle such 
that the suction surface and pressure surface 
velocity distributions close inside the trailing 
edge at perhaps 85 to 90 percent of chord, was one 
possible means for accounting for the effect! of it 
rounded trailing edge and boundary layer sep< ration 
over the rear portion of the suction surface. For 
controlled diffusion blading, the object is to have 
no boundary layer separation. If this is accomp- 
lished, the deviation angle would be expected to be 
small, and there would be justification for allowing 
the suction surface and pressure surface velocities 
to close at the trailing edge, rather than closing 
earlier. 

The constraint Is defined as a non-dimensional 
difference between velocities on the suction and 
pressure surfaces at the trailing edge mesh line, 
and Is expressed: 


The denominator, 15. 24, was chosen to scale the 
constraint to about order one. Upper bound was set 
at zero and lower bound at -1.2F, The velocity 
difference at the trailing edge could vary between 
an upper bound of zero and a lower bcund of 
-19.0 m/sec, thus permitting some closing Inside 
the trailing edge, 

MODIFICATIONS OF ANALYSIS PROGRAMS 

An Important requirement of the optimization 
method is that the objective and constraint func- 
tions bir continuous functions of the design vari- 
ables. This necessitated certain modifications of 
the analysis programs. 

Modification of Inviscid Code 

Experience with the TSONIC code has shown that 
calculations in the trailing edge region can be 
quite sensitive for some configurations. Orienta- 
tion of the blade, trailing edge radius, and grid 
Intersection points can affect surface velocity 
calculation? at or near the trailing edge station, 
sometimes resulting in spurious behavior. Inaccu- 
rate trailing edge velocities will produce incorrect 
gradients of the trailing edge constraint described 
above, and possibly give misleading violations of 
that constraint, 

The means used to avoid or reduce this tendency 
is to incorporate a mass injection model at the 
trailing edge (15). In this model, tangents are 
formed at the intersection of the trailing edge 
circle with the blade surface, and extended to the 
vertical grid line which forms a tangent with the 
trailing edge circle (Fig. 5). The "wake" Is then 
extended downstream with an orientation determined 
by the downstream whirl boundary condition, Experi- 


ence has shown this modeling to reduce the sensi- 
tivity of the surface velocity calculations in the 
trailing edge region. 

Modifications of Boundary Layer Code 

' There Is presently no agreement concerning tlui 
initial *tate of a boundary layer on a compressor 
stator blade in the real flow environment. Some 
observers have measured laminar boundary layers, 
while others contend that due to high inlet turbu- 
lence and unsteady effects, a laminar boundary layer 
cannot persist. For the purposes of this study, the 
Question is somewhat academic. An optimization 
design process can be developed for either case. In 
the present'work the existence of a laminar boundary 
layer Is assumed, which poses the more difficult 
optimization problem. 

The location o f laminar separation and turbu- 
lent reattachment is of crucial importance to the 
optimization search process. The suction surface 
velocity distribution provldeo as input to the 
boundary layer calculation might ideally appear as 
represented in Fig, 6(a). In reality it might 
appear as in Figs. 6(b) and (c), due to the inter- 
relationship of geometric variables such as blade 
stagger, solidity, camber distribution, thickness 
distribution, transition location, and maximum 
thickness location. Boundary layer calculations are 
inltiatoo with a laminar boundary layer, which would 
usually persist to point A. Laminar separation, 
rather than normal transition, occurs there in all 
cases because of the steep adverse pressure gradi- 
ent, Conservation of momentum is assumed through 
the laminar reparation region, with the turbulent 
boundary layer reattached at the next calculating 
station, Turbulent separation is assumed to occur 
when the incompressible form factor reached a criti- 
cal value. 

As originally modelled, point A (Fig. 6) is 
identified as the station at which skin friction 
becomes negative. Any sensitivity to design vari- 
ables can cause a discontinuous jump in point A 
location. This effect carries through to directly 
influence turbulent boundary layer separation loca- 
tion and the objective function. To establish a 
consistent and Conservative criterion, the following 
procedure was coded. Using Lagrangian inter- 
polation, three additional points are placed between 
each station in the high gradient region of the 
velocity vs. distance array. A search procedure is 
begun from the trailing edge region, and locates the 
maximum velocity at the beginning of the high 
gradient region, point B in Fig. 6. Laminar 
separation and turbulent reattachment is effected at 
point 6, 

In addition to the modifications dismissed 
above, several modifications were required relating 
to turbulent boundary layer separation, A separa- 
tion criterion common to compressor blade analyses 
which use integral boundary layer methods is the 
incompressible form factor, Hi. Values of 1.8 to 
2.6 have been proposed and used in the past (e.g., 
von Ooenhoff and Tetervin, Ref. 16). A value of 2.0 
is somewhat conservative and, in the experience of 
the author, has proven to be useful. The program 
was modified to use 2.0 as the critical incompressi- 
ble form factor. 

In normal operation, when a form factor at a 
given station exceeds 2.0, separation is assumed to 
have occurred at that station. If calculation 
stations are 5 percent of chord opart, separation 
location becomes a discontinuous function, changing 
with distance in 5 percent jumps. To correct this. 
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linear Interpolation is used between stations to 
obtain the percent chord location corresponding to 
H< - 2.0. 

Because of the relation between blade angle 
distribution and thickness distribution, the in- 
compressible form factor quite often resembles 
Fig, 7. A maximum form factor can be observed at 
C. A more conventional form factor distribution Is 
also depicted In Fig. 7, where the maximum value Is 
Identified as 0. A search procedure was added to 
locate the maximum form factor, formax, which Is one 
term in the objective function. 

It was observed that allowing turbulent re- 
attachment at a momentum thickness equivalent to 
momentum thickness at laminar separation' often 
resulted in Initial turbulent momentum thickness 
Reynolds numbers less than 320, the minimum value 
experimentally observed for a turbulent boundary 
layer (17). Therefore, as a final modification, 
the code was altered to provide a minimum thickness 
equivalent to a Reynolds number of 320. 

DESIGN RESULTS 

The optimization history Is shown In Fig. 8, 
Most Improvement occurred In the first two Itera- 
tions. At the end of two Iterations a blade had 
been found with no boundary layer separation 
(XSEPOX - 1.0), Reduction of the objective func- 
tion for subsequent iterations Involved reduction of 
FORMAX only, since XSEPOX remained 1.0. All Im- 
provement beyond iteration 2 provided more safety 
margin from the theoretical separation condition, 

CPU time on an IBM 370/3033 for the eight iterations 
was 49.48 minutes. A total of 85 calls on the 
analysis programs were made. 

The initial and final blade shapes, surface 
velocities, and suction surface boundary layer form 
factor are presented and compared in Figs. 9 to 11. 
The pressure surface boundary layer form factor Is 
presented In Fig. 12. 

In the course of optimization, the geometric 
transition location moved forward from 27.3 percent 
of chord to 24.1, and the level of KTC (blade angle 
at transition) shifted downward from 36.7 degrees to 
34.8 (Fig. 9). The maximum thickness location moved 
rearward from 48.2 percent of chord to 53,6. AH 
polynomial coefficients describing the blade angle 
distribution were altered, as would be expected, 
since they were design variables. The polynomial 
coefficients describing the thickness distribution 
were not altered, since they were not design vari- 
ables. However, since the maximum thickness loca- 
tion Itself changed, the actual distribution of 
thickness was altered, as is evident from Fig. 

9(b). If difficulties in achieving a satisfactory 
design had neen experienced, the polynomial coeffi- 
cients for thickness distribution could have been 
addeo as additional design variables, but at the 
cost of increased computing time. Outlet blade 
angle, KODR, changed little during the process. 

Large excursions in KOCR were prevented because It 
is closely related to the trailing edge constraint 
(Eq. (6)). 

The changes effected in the surface velocities 
by the optimization procedure (Fig. 10) are a bit 
more dramatic in appearance than are the geometry 
changes. The peak velocity on the suction surface 
was reduced, as was the large velocity diffusion 
over the front portion of the pressure surface. The 
unconventional waviness of the pressure surface 
velocity is due to the aft location of the maximum 
thickness. Fitting the thickness distribution 


through this maximum thickness location, In combina- 
tion with the forward transition location, results 
In a region of reversed curvature on the pressure 
surface near the maximum thickness location, and Is 
evident on Fig, 9(c). The effect on the flow 
carries across the channel and appears as a small 
wave on the suction surface as well. Aside from the 
dubious aesthetic appearance, no adverse aerodynamic 
effects can be attributed to this behavior. The 
calculated boundary layers appear well-behaved, with 
the maximum Incompressible form factor on the suc- 
tion surface being 1.924, and on the pressure sur- 
face 1,780. 

In completing the design of the stator, only 
one other blade section was optimized, the hub sec- 
tion at the Inner endwall. This blade element will 
be In the wall boundary layer, so that true two- 
dimensional flow Is not expected to exist. Result- 
ing transition location and maximum thickness loca- 
tion wore not greatly different from the values 
found at the 90 percent span location. Blade angle 
polynomial coefficients different from those ob- 
tained at 30 percent span were obtained. However, 
for reasons relating to the blade stacking proce- 
dure, which will be described below, the polynomial 
coefficients obtained at the 90 percent span loca- 
tion were used also at the 100 percent span. The 
resulting two-dimensional calculations for the blade 
with these coefficients Indicated no boundary layer 
separation. 

All other blade sections were specified based 
on the optimized design obtained at the 90 percent 
span section. Each of these blade sections, which 
He on streamlines, 1 j then radially stacked. 
Fabrication coordinates are interpolated at several 
planes parallel to the axis of rotation of the com- 
pressor. In principle, each of the blade sections 
on the six chosen streamlines utilized could be 
designed by optimization. This could and probably 
would result In six different sets of transition 
location, maximum thickness location, and blade 
angle polynomial coefficients. The fabrication 
coordinates are generated by a design point stream- 
line curvature code. As Input to this code, a 
radial curveflt of each of the polynomial coeffi- 
cients must be provided. Transition location and 
maximum thickness location for each blade section 
are input directly. The blade coordinates on each 
streamline section are then generated. Finally, 
coordinates at the horizontal fabrication planes are 
obtained by interpolation, based on a cubic fit of 
the blade coordinates at the four stri amlines most 
closely straddling the desired fabrication plane. 

Because of the curvefitting at various stages 
of this process, prudence suggests avoiding the 
possibility of large radii! variations in the design 
parameters. Therefore, a constant radial distribu- 
tion of each parameter was sought, with one excep- 
tion. Maximum thickness location was arbitrarily 
moved forward to 47 percent of chord for all sec- 
tions between the tip and 70 percent of span from 
tip. Although it was not necessary to do this, the 
effect was to relieve the reversed curvature condi- 
tion on the pressure surface. The transition loca- 
tion and all polynomial coefficients were maintained 
at the same values obtained for the optimized 90 
percent span section. At the 100 percent span sec- 
tion, transition and maximum thickness locations 
found from optimization at that section were used 
(0.26 and 0.52 respectively), and blade angle poly- 
nomial coefficients equivalent to those at 90 per- 
cent span were used. Thus, neither the polynomial 
coefficients for blade angle nor maximum thickness 


varied radially. Transition location was constant 
from tip to 90 percent span, and differed only 
slightly at 100 percent span. Maximum thickness 
location was constant from the tip to 70 percen', 
span at 0.47, moved rearward to 0.53 at 90 percent 
span, and slightly forward to 0.52 at 100 percent 
span. The exit blade angle, KOCR, varied only 
slightly from tip to hub In a range from -3.7 to 
-4.0. Since design exit flow angle Is zero degrees 
for all sections, the negative value of KOCR re- 
presents deviation angle. If, Indeed, the boundary 
layer does not separate from the blade as theoreti- 
cally predicted, the deviation angles of about 
4 degrees may be more realistic than they appear to 
be, The blade geometry and surface velocity dis- 
tributions for the blade sections at midspan and 10 
percent span from tip are shown In Figs. 13 and 14. 

SUMMARY AND CONCLUDING REMARKS 

A method has been presented for automated com- 
pressor blade design using numerical optimization 
technique*. The method was applied to the design of 
a control ^-diffusion stator blade row. T>-' ee 
analysis programs were coupled to the numerical 
optimisation program: a blade geometry generation 
program which uses polynomial representation for 
blade angle and thickness distributions, a com- 
pressible, Invlscld flow program, and an Integral 
boundary layer program. Seven of the nine design 
variables were related to blade angle distribution, 
another located the maximum thickness of the blade, 
and the last controlled camber and deviation 
angles. Two constraint functions operated In the 
geometry program to produce shapes with controlled 
diffusion velocity olstrlbuttons. Constraint func- 
tions applied In the flow analysis programs limited 
suction surface velocities to the subsonic regime, 
limited the velocity diffusion on the pressure sur- 
face, and set the trailing edge condition for 
Invlscld calculations. The objective function, 
which was minimized, was of a penalty function form, 
and effectively produced a blade whose suction sur- 
face turbulent boundary layer did not separate. 

The optimization procedure for the subject 
blade section required eight major Iterations in- 
volving 85 calls on the geometry/ aerodynamic analy- 
sis programs. Total CPU time on an IBM 370/3033 
computer was 49. <*5 minutes. 

When using the numerical optimization proce- 
dure, It was essential that the gradients of the 
constraint and objective functions be smooth and 
accurate. Therefore, some modifications of the 
analysis programs were necessary to ensure that 
these functions were continuous. 

The design problem, as formulated here, pro- 
duced a blade shape which satisfied the design 
criteria, while holding the polynomial coefficients 
describing thickness distribution, the value of 
maximum thickness., and the incidence angle con- 
stant. The method thus still offers great flexi- 
bility for adaptation to more demanding design 
requirements. 
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(c) Blade shape. 

Rgure 2. - Initial blade geometry used In optimization search. 




Figure 3, * Initial surface velocity distribution. 



Figure 4. - Initial design-suction surface turbulent 
boundary layer form factor distribtion. 
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Figure 7. - Turbulent Incompressible form factors for different classes 
of blades. 
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Figure 10. - Comparison of initial and final blade sur- 
face velocity distributions. 






